Is a "homogeneous" description of dynamic heterogeneities possible? 
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We study the simplest model of dynamic heterogeneities in glass forming liquids: one-spin fa- 
cilitated kinetic Ising model introduced by Fredrickson and Andersen [G.H. Fredrickson and H.C. 
Andersen, Phys. Rev. Lett. 53, 1244 (1984); J. Chem. Phys. 83, 5822 (1985)]. We show that the 
low-temperature, long-time behavior of the density autocorrelation function predicted by a scaling 
approach can be obtained from a self-consistent mode-coupling-like approximation. 

PACS numbers: 64.70.Pf, 64.60.Ht, 64.60.Cn, 05.50+q 



Recently, a new general paradigm for dynamics of su- 
percooled liquids has been proposed^. It focuses directly 
on the spatial and temporal structure of dynamic het- 
erogeneities, i.e. regions of higher-than-average mobility. 
The essence of this new approach is a mapping of a mi- 
croscopic model of a supercooled liquid onto an effective 
model describing dynamics of the heterogeneities^. It has 
been argued on physical grounds that very simple effec- 
tive models can capture the most important features of 
the heterogeneities' dynamics. In the initial worki two 
models were used: a one-dimensional one-spin facilitated 
kinetic Ising model, i.e. the Fredrickson- Andersen (FA) 
model 3 , and its asymmetric version, i.e. the East model 4 . 
In the subsequent work£, more complicated (and more 
realistic) models were used. However, it is the simplest 
model (i.e the FA model) that is most often used-^ to 
illustrate the main points of the new paradigm and to 
contrast it with other approaches to the glass transition 
problem like the mode-coupling theorjs& and the land- 
scape paradigm 9 . 

The simple models describing dynamic heterogeneities 
have trivial static properties. However, they exhibit in- 
teresting dynamics. Hence, there arises a new prob- 
lem: how to study the dynamics of the effective mod- 
els. One should notice that solving this problem will 
be, most likely, considerably simpler than solving the 
original, microscopic model of a supercooled liquid: the 
effective models used so far are lattice-based and thus 
have a greatly reduced number of degrees of freedom; 
even more importantly, these models usually involve a 
small parameter (i.e. density of heterogeneities (mobile 
regions)) which can facilitate theoretical analysis. 

The original work 1 concentrated on the analysis of the 
structure of trajectory spaces of the FA and East models; 
it was argued that the glass transition coincides with an 
entropy crisis in trajectory space rather than in config- 
uration space. This is an ingenious suggestion because 
it focuses directly on the system's dynamics and thus 
it should be applicable to other glass-like systems with 
trivial equilibrium propertiesifliii. 

In Ref. 1 , in order to get explicit results for dynamics 
of the FA and East models, the random walk theory was 
combined with scaling arguments. Subsequently 2 a dy- 
namic field-theoretical model has been derived and ana- 
lyzed using the renormalization group (RG) methods. It 



was argued that the d-dimensional FA model belongs to 
the same universality class as the field-theoretical model. 
RG analysis showed that, for d > 2, there is no finite 
temperature phase transition and the low temperature 
scaling properties are determined by the directed perco- 
lation critical point. 

In our opinion, the simplicity of the effective mod- 
els suggests that it might be possible to analyze their 
dynamics using well-known methods of non-equilibrium 
statistical mechanics. Such an approach was utilized by 
Fredrickson and Andersen^: the time-dependent density 
correlation function was expressed in terms of a mem- 
ory function; the latter function was then analyzed us- 
ing an exact diagrammatic expansion 1 ^. The resulting 
self-consistent equation for the density correlation func- 
tion predicted a power law dependence of the relaxation 
time on the density of heterogeneities. However, neither 
Ref£ nor later work on the FA modeli 3 . reproduced cor- 
rect scaling exponent for the relaxation time. This was, 
perhaps, the reason for the recent claimi that no "ho- 
mogeneous" , mean-field— approach is able to reproduce 
scaling predictions for the FA model. 

In this note we present a simple theory that does re- 
produce previously predicted scaling behavior— of the 
FA model density correlation function. The theory re- 
lies upon a self-consistent factorization approximation of 
a pair-density correlation function into a product of two 
density correlations. This approximation is analogous 
to one used in the mode-coupling theory^; the impor- 
tant difference is that we use the factorization approxi- 
mation for pair-density involving non-nearest-neighbors. 
We argue that while this makes the structure of the the- 
ory somewhat complicated, it is necessary from a physi- 
cal point of view since the dynamics of nearest-neighbor 
up-spins are drastically different from that of individual 
spins. 

The theory presented here shows that it might be pos- 
sible to use established "homogeneous" methods to ana- 
lyze dynamics of effective models of heterogeneities. 

The FA model consists of a chain of N non-interacting 
Ising spins in a unit magnetic field. We use density vari- 
ables rii (n, = 1 corresponds to zth spin pointing up; 
physically, it corresponds to a more mobile region located 
near site i). In terms of rii the Hamiltonian is H = n,. 
The average density is (nj) = c = l/(l+exp(l/T)) where 
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T is the temperature. 

We consider the normalized time-dependent density 
correlation function, dj(t), 



Cij(t) = (Srii exp (fit) Srij) I ((Srii) 



(1) 



where Srii is the deviation of the density at site i from 
its equilibrium average, Sni = rii — c; note that Cjj (t — 
0) = Si j. In Eq. O is the evolution operator; time 
evolution consists of single spin-flip dynamics with rate 
that is proportional to the number of nearest neighbors 
in the spin up state. Explicitly, 



n 



N 
i=l 



(1 - Si)u>i 



(2) 



where Si is the spin- flip operator, Sifii = — rii + 1, m, is 
the number of nearest neighbor up-spins, m, = Rj+i + 
rii_ i, and Wi is the factor maintaining detailed balance, 
Wi = c + n,i(l — 2c). 

Some earlier simulational studies^ and recent com- 
putational and theoretical worksi^ used the so-called 
persistence function to characterize dynamics of the FA 
model. The persistence function is defined as a fraction 
of spins that remain unflipped after a given time. While 
this function is easier to determine from simulations than 
the density correlation function, it is somewhat awkward 
to analyze theoretically because it is defined via a time- 
nonlocal constraint. 

The starting point of our analysis is the memory func- 
tion representation of the density correlations function, 

£ fa + Mlj r (z) (dn^Sn,)- 1 ) (zC jtk (z) - 6 j>k ) 



= (SriiftSni) (SniSni) 1 C ijfc (z) 



(3) 



Here Ci.j (z) is the Laplace transform of Cij (t), Ci.j (z) = 
J* °° dt cxp(— zt)dj(t), and M}j r (z) is the Laplace 
transform of the irreducibleA£ii& memory function. The 
memory function can be expressed in terms of the pair- 
density correlation function involving nearest neighbor 
spins, 

Mly(t) = Ai d (t) + A^ ltj (t) + Aij-xQt) + Ai.ij.^t). 

Note that if the M}j r is neglected, the density correla- 
tion function Cjj is diagonal; motion of the isolated up- 
spins is possible only due to the coupling to pair-densities. 

The pair density correlation function Aij (t) evolves 
with 1-spin irreducible evolution operator f2 llrr , 



Aij(t) = (SriiSrii+i exp (S} Urr t) SnjSnj+i) (5) 



where 



+ niiSrii) (SniUiiSni) 1 (Sniirii 



Qv 
(6) 



In Eq. JBJ Qi = 1 — P\, and P% is a projection operator 
on the one-spin density subspace, 



A = 5n i) (( S 



{Srii 



(7) 



It should be noted at this point that the standard se- 
quence of approximations leading to a mode-coupling 
theory- would start from Eq. JSJ. Briefly, the pair- 
density correlation function would be factorized into a 
product of two density correlations, and the irreducible 
operator Q ltrr would be replaced by niS. However, a 
closer investigation shows that the time dependence of 
Aij is very different from that of the product of CijS. 
In particular, if in an equation of motion for Aij one 
neglects the memory function term, one gets a (lattice) 
diffusion equation. This is in contrast to Eq. J3J that, 
upon neglecting the memory function term, leads to a 
site-diagonal density correlation function. 

The strategy followed in this note is different: we pro- 
pose to continue using the memory function approach 
until we generate a function that can be expressed in 
terms of pair-densities involving non- nearest neighbors, 
in terms of e.g. SriiSrii+i, I > 2. With this goal in mind, 
we derive a memory function representation for the pair- 
density correlation function, Aij. It turns out that the 
2-spin irreducible memory function that enters the equa- 
tion of motion for A^j can be expressed in terms of the 
triple-density correlation function, Pij, 

Mfj r (t) = B itj (t) + Bi_ij(t) + Bij-i(t) + Bi-u-iW, 

(8) 

where Bi^ is defined as 

Bij(t) = (5niSni+iSrii+2 exp {fl 2lrr t) SnjSnj + iSnj + 2) , 

(9) 

and the 2nd order irreducible evolution operator f2 2lrr 
is given by a formula analogous to JSJ- It should be 
noted that £l 2lrr involves operators Q2 projecting on the 
space orthogonal to the density and to the pair-density 
involving nearest neighbor sites, Q% = 1 — Pi — P%. Here 
Pi is given by Eq. {7J and P2 is defined as 



Pi = ^2 $ n id n i+i) ({SniSn i+ if\ (SriiSn 



(10) 



Next, we derive the memory function representation 
for Bij. Its irreducible memory function is given by the 
following expression 



M-J r {t) = (SnMi+iSni+2 (m i+2 + m.) Q3 

exp (fl 3irr t) Q 3 (to,-+2 + n%j) Su ,r)n , . Su , . 2 ) , (11) 



where Q3 = 1 — Pi — P% — P3, and P3 is defined as 



P 3 = SniSm + iSn i+ 2) ((SnjSrii+iSni+zy 

i 

(SniSni + i5n i+ 2. (12) 
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Finally, the 3rd order irreducible evolution operator tt 3lrr 
in Eq. Ijlll) is given by a formula analogous to iJBJ. 

At this point we note that Mf"' r is a correlation func- 
tion of a quantity that has a non-zero overlap with the 
pair-density involving non-nearest neighbor sites: 

i5ni5n i+2 Qz (m i+2 + m,) Sn^n.^Sni^ 

= 2(c(l-c)) 3 . (13) 

Thus, it is at this stage that we propose to follow the 
usual sequence of approximations 8 . As a result we get 
the following expression for the memory function M 3 J r 



in terms of the density correlations Cjjs: 

M 3 J r (t) « 4 (c(l - c)) 4 C id (t)C i+2J+2 (t). (14) 

The final self-consistent equation for the density cor- 
relation function is easiest to present in terms of the 
Laplace-Fourier transform, C(q; z), 

C(q;z)=J2^M-iQk)C k , (z). (15) 

It reads 



C(q-z) 



z + 



2c 

2(1 -c) cos 2 (<?/2) 



(16) 



z + 2(l - c) sin 2 (g/2) + 



2c 



2(1 - c) cos 2 (<?/2) 



z + 2(l-c) sin 2 (g/2) + 



l + 2c(l-c)<M<?;z) 



where (f>2(q;z) is the Fourier-Laplace transform of 

Ci j (t)Ci+2j+2 (t) ■ 

To investigate the low temperature (i.e. low c) behav- 
ior of the density correlation function we consider the 
scaling function C(q*;z*), 

C(q*;z*) = lim c 3 C(q*c; z*c 3 ) 



dt*exp(-z*t*)C(q*;t*) (17) 



where C(q*;t*) = lim c _>o C(q*c\ t*c~ 3 ). It can be easily 
shown that in the limit c — > the self-consistent equation 
Hlfijl reduces to the following c-independent equation for 
the scaling function C(q*; z*): 



C(q*;z*) = 



(18) 



z* + q* 2 /2 + l/Mq*;z*y 
Here (f>2(q*;z*) is the Laplace transform of 4>2(q*]t*), 
(j>2(q*;t*) = lim c~ 1 (j)2{q*c; t*c~ 3 ) 

?—C(k*;t*)C(q*-k*;t*) (19) 

Z7T 

The scaling (II Til and c-independence of Eq. (|18|l shows 
that, in the low c limit, up-spins (i.e. mobile regions) 
undergo diffusive motion with a diffusion constant D oc 
c; up-spin diffusion is cut-off at long times by a longest 
relaxation time r oc c~ 3 . These results agree with known 
low c scaling predictions^. 



In order to get quantitative results for the diffusion 
constant and the relaxation time one would have to solve 
the scaling equation l|18[) combined with Eq. I|19fl . The 
numerical solution of Eqs. I|18ll9fl an d a comparison of 
the resulting density correlation function with computer 
simulations are left for a future study. Here we note that 
a single iteration of Eqs. (|18I19|) starting from a zeroth 



order approximation &°\q*; z" 
leads to 

C (1) (q*;z*) 



) = 1/ (z* +aq* 2 /2 + < 



z* + q* 2 /2 + (Aaz* + 8ab + a 2 q* 2 ) 1/2 ' 

(20) 

Whereas the zeroth order approximation has a sim- 
ple pole, the first order approximation has a pole 
and a branch cut. The simplest self-consistency condi- 
tion that allows one to estimate a and b (and thus the 
diffusion constant and the longest relaxation time, re- 
spectively) is to identify the pole in and the one in 
C*W. This leads to the following approximate results for 
the diffusion constant^ and the longest relaxation time: 
D = 0.4c, t = 0.3125c- 3 . 

One should note that the structure of Eq. l[T%|l is sim- 
ilar to that of the self-consistent equation of the mode- 
coupling theory of the glass transition^. The very impor- 
tant difference is the q* 2 /2 term in the denominator. It is 
this term that cuts off the ergodicity-breaking transition 
that is generic for mode-coupling theories. 

In summary, we have shown here that a simple self- 
consistent, mode-coupling- like approximation reproduces 
known scaling predictions for the density correlation 
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function of the FA model. It would be of interest to 
investigate whether a similar approach can be applied 
to more complicated models like the East model or the 
Kob- Andersen lattice gas model 21 , which recently have 
attracted renewed interes1*i*2£. 



Note added: After this paper had been submitted 
a preprint presenting an RG approach to the one- 
dimensional FA model appeared^. 

Support by NSF Grant No. CHE-0111152 is gratefully 
acknowledged. 
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